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Abstract 

Zero-bin subtractions are required to avoid double counting soft contributions in collinear loop 
integrals in Soft-Collinear Effective Theory (SCET). In traditional approaches to factorization, 
f-^. \ double counting is avoided by dividing jet functions by matrix elements of soft Wilson lines. In 

this paper, we compare the two approaches to double counting, studying the quark form factor and 
■ deep inelastic scattering (DIS) as xb — ► 1 as examples. We explain how the zero-bin subtractions 

in SCET are required to reproduce the well-established factorization theorem for DIS as xb — > 1. 
1} . We study one-loop virtual contributions to the quark form factor and real gluon emission diagrams 

in DIS. The two approaches to double counting are equivalent if dimensional regularization (DR) 
is used to regulate infrared (IR) divergences. We discuss in detail ambiguities in the calculation 
of one-loop scaleless integrals in DR in SCET and perturbative QCD. We also demonstrate a 
nontrivial check of the equivalence of the zero-bin subtraction and the soft Wilson line subtraction 
{Sj . in the virtual two-loop Abelian contributions to the quark form factor. 
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I. INTRODUCTION 



High-energy processes with hadrons in either the initial or final state are amenable to 
perturbative QCD (pQCD) treatment as long as there is a hard scale, Q, much larger than 
Aqcd- The essential idea of pQCD factorization [l| is to separate the different relevant 
scales into a well-defined quantities that capture the physics at the given scale. By doing 
so one is then able to distinguish the perturbative short-distance contribution from the 
nonperturbative long-distance one. 

In order to properly establish such factorization of scales it is clear that one has to avoid 
double counting among the functions that appear in the factorized physical quantity. In gen- 
eral these nonperturbative functions receive leading contributions (in powers of Aq CD /Q 2 ) 
from specific regions in momentum space. For example, the contribution to a jet function 
may be dominated by partons with momenta whose components scale as 



where A characterizes the off-shellness of the partons in the jet function, and A 2 ~ QAqcd or 
Aq CD , depending on the process considered. In principle, one could restrict all momenta ap- 
pearing in loop calculations to satisfy the scaling in Eq. ([T]) , but this is impractical. Instead 
one typically integrates over all momentum space in loop integrations, making approxima- 
tions appropriate for the assumed scaling. This practice is well-known in, for example, the 
method of regions [2], where approximations appropriate to a specific momentum region are 
performed at the level of the integrand, but the integrals are over all momentum space. In 
cases where contributions from other momentum regions are not sub-leading, one encounters 
double counting. 

Consider, for example, deep inelastic scattering (DIS) in the end point region, xb — ► 1, 
where xb is the well-known Bjorken variable. In this region, Sterman [3[ has shown that 
to all orders in perturbation theory one can (re-)factorize the non-singlet DIS structure 
function in moment space according to the formulae 



Here F 2 ^{Q 2 ) is the Nth moment of the structure function and xb — ► 1 implies N ^> 1. H 
is the hard contribution, 0at is the Nth moment of a modified parton distribution function 
(PDF), Jjv is the Nth moment of an outgoing jet function, and Sn is the Nth moment of a 
soft factor representing the emission of soft gluons. Explicit and gauge invariant expressions 
for the PDF, jet, and soft functions will be given below. In the factorization theorem of 
Ref. [§], the jet, soft, and parton distribution functions are not defined in a gauge invariant 
way, and through a judicious choice of gauge one can eliminate double counting. However 
when one defines the same quantities using Wilson lines to make them gauge invariant, 
the factorization theorem breaks down and has to be modified to take into account the 
double counting of soft momentum modes in the soft factor, Sn, and in the collinear matrix 
elements, and <p^. The Wilson lines in the matrix elements are the cause of the double 
counting problems. In Ref. it was proposed that for gauge invariant quantities the 
factorization theorem should take the form 
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A similar factorization theorem with a soft factor subtracted from the transverse momentum 
dependent PDF and jet function was also proposed by Ji, et. al. [5] for semi-inclusive DIS. 
Ref. 0] checked the theorem explicitly at one-loop and gave arguments for the theorem to 
hold to all orders in perturbation theory. Another well studied quantity in pQCD is the 
quark form factor and its factorization into jets and soft contributions. The double counting 
problem for this quantity is studied in Refs. @, @], where it is argued that the same soft 
factor has to be subtracted from each one of the two collinear jets, in a manner similar to 
Eq. (j3J), to avoid double counting. 

Recently, the problem of factorization has been addressed in an effective field theory 
approach using Soft Collinear Effective Theory (SCET) [1,0]. In this approach to factor- 
ization, QCD is matched onto SCET, a theory containing collinear and soft modes whose 
momenta scale as 1 

collinear ~ Q(1,A 2 ,A) 

soft ~ Q(A 2 ,A 2 ,A 2 ). (4) 

After matching QCD onto SCET, a field redefinition [gj can be used to decouple the soft and 
collinear modes at the level of the Lagrangian. In SCET, jet functions are matrix elements 
of collinear fields and soft functions are matrix elements of soft fields. Note that in this 
paper the terms soft and collinear are defined by Eq. pi) . For example, if a gluon whose 
momentum components all scale as 0(QX 2 ) is emitted from a jet with small angle we still 
refer to it as a soft gluon. 

The subject of double counting has gained renewed interest in the context of SCET due to 
the work of Manohar and Stewart [13] • In SCET, power counting is made manifest using the 
label formalism [ll[. Let <f>(x) denote a generic full theory field, the corresponding collinear 
field in SCET, 4>p{ x )i is defined by 

= 5>-**&(a0 . (5) 

The p are 0(Q,QX) label momenta which correspond to the large parts of the collinear 
momenta, while derivatives on (f>p{x) give the 0(QX 2 ) residual momenta. Loop integrations 
involve both a sum over labels and an integral over the residual momentum which in practical 
calculations are combined into an integral over all momentum space, 




But this raises the issue of double counting, as the integral over all momentum space in- 
cludes the region in which the label momentum vanishes - the "zero-bin" - which is already 
included in the soft sector of the theory. To correctly calculate the collinear graphs, one 
must subtract the zero-bin contributions. The zero-bin subtraction is necessary for the 
proper interpretation of 1/e poles in SCET loops [13], and also necessary so that collinear 



In Ref. Um, 

modes with all momentum components 0(QX) are called soft, while modes with all momen- 
tum components 0(QX 2 ) are called ultra-soft, or usoft. The 0(Q\) modes play no role in this paper, so 
we will neglect them, and we will refer to 0(Q\ 2 ) modes as soft rather than ultra-soft. 
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plus soft real gluon emission graphs correctly reproduce real gluon emission QCD diagrams 



in the appropriate regions of phase space [10|, ll2| 



The zero-bin subtraction in SCET is clearly related to the soft subtraction in Eq. ([3]). 
In calculating the zero-bin contribution to a particular Feynman graph, one must apply 



SCET power counting. In the minimal zero-bin subtraction scheme of Ref. [10| the zero- 
bin contribution is only removed if its contribution scales as A , i.e., is leading order in 
SCET power counting. Only the zero-bin of the 0(Q\ 2 ) component of the gluon couples to 
collinear fields at 0(A ) so in this scheme collinear couplings to other zero-bin modes can 



be neglected. Exploiting this fact, Ref. [13| shows that the difference between calculating 



collinear matrix elements naively and including the zero-bin subtraction amounts to using 
two different collinear Lagrangians that are related by a field redefinition of the collinear 



fields [13]. This field redefinition is similar to the field redefinition that is used to decouple 
soft modes from collinear modes to prove factorization in SCET. The field redefinition can 
be used to decouple the zero-bin mode from collinear modes at the price of inserting zero-bin 
Wilson lines into the collinear matrix elements. Thus, the naively evaluated collinear matrix 
element factorizes into a properly evaluated collinear matrix and a matrix element of zero- 
bin Wilson lines, establishing the equivalence of zero-bin subtractions and soft Wilson line 
subtractions at lowest order in A. At higher orders in A, zero-bin modes of other components 
of the gluon field are important so it seems that the equivalence only holds to lowest order 
in A. 

The purpose of this paper is to study the equivalence between zero-bin subtractions and 
the soft Wilson line subtractions in more detail. We focus on the on-shell quark form factor 
and the factorization theorem for DIS as xb — » 1- One important issue is the choice of 
IR regulator in the calculation of the zero-bin contribution. The field redefinition which 
decouples the soft or zero-bin modes and collinear modes will only leave on-shell matrix 
elements invariant, and therefore off-shellness is not a suitable regulator [13]. Soft Wilson line 
subtractions and zero-bin subtractions are not equivalent if offshellness is used to regulate 
the IR divergences. 

In our calculation of the quark form factor, we will use DR to regulate both ultraviolet 
(UV) and IR divergences. The individual one-loop soft diagrams, the collinear diagrams and 
their zero-bin subtractions are all ill-defined, and only their sum yields a well-defined result 
which reproduces the IR divergences of QCD. Fortunately, in this case the equivalence of 
the soft and the zero-bin subtractions can be easily established at the level of the integrands. 
At higher orders the zero-bin subtraction involves an iterative procedure which requires the 
subtraction of graphs in which some lines have collinear momentum and some have soft 
momentum, so demonstrating the equivalence of soft and zero-bin subtractions at higher 
orders is somewhat subtle. For example, one has to be careful about the order in which 
limits are taken to correctly compute the zero-bin subtraction at two loops. Another issue is 
that the equivalence of zero-bin and soft subtractions does not hold at the level of individual 
Feynman diagrams, but only for the sum of all diagrams at a given order. In some cases 
this requires a cancellation between the zero-bins that appear in different diagrams. In this 
paper, we explore these issues by calculating the zero-bin subtractions in the two-loop abelian 
diagrams contributing to the quark form factor. We also show that they are equivalent to 
soft subtraction. 

Another important point of this paper is to emphasize the importance of the zero-bin 
subtractions in the collinear matrix elements that appear in the SCET derivation of the 
endpoint factorization in DIS as xb —>■ 1. We will assume a factorization theorem of the 
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form of Eq. (j2J) where the definitions of Sn, 4>n, and Jn in terms of SCET fields will be 
given below. We will show that the one- loop zero-bin contributions to <$n and Jn are equal 
to Sn, which confirms Eq. We will also show that at one-loop Eq. reproduces the 
large N limit of the full QCD calculation of DIS as Xb — > 1- In these calculations we again 



use DR to regulate the IR divergences. Our results are consistent with those in Ref. [15 
however, in our treatment all Wilson lines are defined on the light-cone. For other SCET 
treatments of DIS in the xb — » 1 limit, see Refs. 0,HE3- 



This paper is organized as follows. In section II, we discuss the space-like quark form 
factor in SCET and present one-loop results for the various contributions. Once zero-bin 
contributions are included it is straightforward to verify that SCET reproduces the IR 
divergences of QCD. In section III, we give our analysis of DIS at x — > 1 in SCET. In 
section IV, we study the zero-bin subtractions for two-loop virtual abelian contributions to 
the quark form factor and demonstrate agreement with the soft Wilson line subtractions. 
In section V we conclude. 



II. ZERO-BIN AT ONE-LOOP: QUARK FORM FACTOR 

In this section we calculate the quark form factor to one loop using DR to regulate both 
the UV and IR. It is straightforward to see that the zero-bin subtractions are reproduced by 
the soft subtractions for this choice of regulator. The quark form factor may be the simplest 
quantity to analyze in SCET, however it serves well to establish some of the relevant issues. 
We take the incoming quark to be moving along the +z direction so its four-momentum is 
Pi = (PiiPi iPi) = = 73 n/1 - ^he i ncom i n g quark scatters off a space-like photon 

whose invariant mass is q 2 = — Q 2 < 0. The final state quark is moving in the —z direction 
so its momentum is pr, = iptiP2 1P2) = (0> ~j|> 0) = 75^- 

The full QCD electromagnetic current, j M = ip(x)'j' J "ip(x), is matched onto an effective 
one, 

f = C(Q 2 )j» ff = C(Q 2 ) InW^W^n 

= c(Q 2 ) S Vfy fl Vy„wJ°)t^ . (7) 

The manipulations performed in Eq. are standard in the SCET formalism. The form of 
the current is fixed by collinear gauge invariance. The second line is a consequence of making 
the field redefinition of the collinear quark and gluon fields that decouples soft gluons from 
the collinear sector of the leading order SCET lagrangian. The collinear Wilson line, W n , 
and the soft Wilson line, Y n , are given in momentum space by 

^— ' ml n ■ qin ■ [qi + q 2 ) ■ ■ ■ n ■ (2^-1 Qi) 

m=0 perms v ' vi — <t—± * / 

and 



Y = y~9) m n ■ A s - ■ - n ■ A s 

n m! n ■ qin ■ (g x + q 2 ) ■ ■ ■ n ■ QXi Qi) ' 

m=0 perms ^ — n—± 



A n ^ q . is a collinear gluon field and all n ■ q$ scale as Q. A s stands for soft gluon field and n ■ q^ 
scales as QX 2 . 
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FIG. 1: The quark form factor in SCET. In the first two graphs the thick line attached to the 
fermion line represents a collinear Wilson line, while in the last graph both thick lines are soft 
Wilson lines. 



The matrix element of j^j- factorizes into 

(q(p2)\f\q(pi)) = C(Q 2 )r[Jout x/x J m x S] , (10) 

where the jet and soft functions are defined by 

J in = (0|Wf )t e (o)| ?(pi)); j out = (qfaftfflwP |0>, S = (0\YiY n \0) . (11) 

In the jet functions the fields are collinear so only fields with nonvanishing labels appear. As 
argued by Lee and Sterman 13) , only the zero-bin of the n ■ A n>q component of the collinear 



gluon couples to ffi at leading order in A. Furthermore, the Wilson line, Wn^\ which 
has no zero-bin gluon is related to one with a zero-bin by multiplication by a Wilson line 
constructed from the zero-bin component of the n ■ A n ^ q gluon. Thus the collinear matrix 
element is equal to the naively evaluated matrix element divided by a matrix element of soft 
Wilson lines: 

(oi^r t ^ ) i g (Pi)> = <0| ^ w|g(Pl)> , (12) 



where 



Y n (x) = Pexp 



o 



ig s I dsn ■ A n fl(ns + x) 



(13) 



A n> Q is the zero-bin of the collinear gluon, and our definition of the covariant derivative is 
£)n _ Q/i _ ig s j[t J '. w e will use the same notation for the soft Wilson line and the zero-bin 
Wilson line because the matrix elements are identical and it is always clear from context 
which is relevant. In Eq. (fT2|) the hats on the fields signify that contains the zero-bin 
gluon and that the Lagrangian includes the coupling of £ n to zero-bin gluons, so this matrix 
element corresponds to the naive evaluation of the collinear matrix element. In our notation 
(0|wi O)t ^ 0) |g(pi)) corresponds to (0| W£f'Jg(pi)) of Ref. 0, {0\W^ n \q( Pl )) corresponds 
to (0|iy""^|g(pi)), Y n corresponds to U^, and Y n corresponds to Q n . 

At one-loop the Feynman diagrams that contribute to the collinear n-jet, collinear n-jet 
and the soft factor are given in Fig. [lj All these diagrams are scaleless in DR and therefore 
all integrals will vanish when one sets £uv = e m- Using SCET Feynman rules to calculate 
the n-collinear diagram and naively integrating over all momentum space, the result is 

I n = -2ig 2 s C F {^Y J (27r)d { _ 2pk - + k 2 + P lQ] { _ k+ + iQ] [k 2 + lQ] « ( 14 ) 
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where p = p\ = p 2 — Note that the collinear graph is not I n but rather I n — I n , where 
I n fl is the zero-bin contribution. We refer to I n as the naive n-collinear contribution. The 
naive n-collinear contribution is similar, 

n = ~ 2tg ° CF ^ > J ( 27r )<* [-2pk+ + k 2 + i0] \-k- + iO] [k 2 + iO] ' { } 

and the soft contribution is 

I s = -ig 2 s C F ^ 2 f f ^ )d{ _ k+ + m ^ +mk2 + iQ) ■ (16) 

Let us now compute the zero-bin subtraction for I n . It is easy to check that the zero-bin 
associated with the virtual fermion line is sub leading in A. The only zero-bin needed is 
associated with the virtual gluon. In this zero-bin region, k + scales as QX 2 instead of Q so 
the numerator is simply p. Furthermore, k 2 scales Q 2 A 4 instead of Q 2 \ 2 thus it should be 
dropped relative to pk~ . Similar arguments hold for I n . The result is simply that 



ln,0 = Infl = -ig 2 s C F (n 2 ) £ J 



d d k 



(2vr) d (-k+ + iO)(-Jfe- + iO)(k 2 + iO) 
= I. (17) 

The complete one-loop SCET contribution to the quark form factor is 

I — {In ~ Infl) + {Ifi — Ififl) + Is 

= I n + I n -I s . (18) 

Note that the net effect of the zero-bin subtractions at one-loop is to reverse the sign of the 
soft contributions. Thus the zero-bin subtraction and soft subtraction of Eq. ([3]) yield the 
same results for the one-loop graphs when they are regulated in DR. 

One subtlety in Eq. fflBl is that the integrals I n , I n and I s are all ill-defined. However, 
in Eq. ([TBI the integrands can be rearranged so that I is expressed in terms of integrals 
which are well-defined, then it is straightforward to show that I reproduces the IR of QCD. 
Eq. (|18p can be written as 



-2ig 2 s C F {^ 



d d k ( 2p[p — k + — At 



{27i) d { \-2pk+ + k 2 + i0][-2pk~ + k 2 + i0][k 2 + iO] 
2 

+ 



[-2pk+ + k 2 + iO] [-2pk~ + k 2 + iO] 

k 2 



2[-2pk+ + k 2 +i0][-2pk~ + k 2 + iO)[-k+ + iO][-k~ + iO] 
h + h + h, (19) 
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where 



h = -2tg 2 C F (^ 



d d k 2p(p- k + - k~) 

(27r) d [-2pk+ + k 2 + i0][-2pk~ + k 2 + i0][k 2 + iO] 



a s ^4W r(i + £lR )r 2 (- £lR ) 

"2^ l^J r(2 - 2, IR ) ' (20) 



J 2 2ig 2 s C F {^) £ J ^ d ^_ 2 _ k+ + k2 + + fc2 + 



- + 2^ {-&-) r (2 - 2euv ) ' (21) 

3 - - iff,OFl/i J y (27r) d 2[-2pA;+ + A; 2 + iO] [-2pfc- + fc 2 + iO] + iO] [-k~ + iO] 

a /VV 1 r(i- guv ) 2 r(i + e uv ) 

- + 2^ Cf 7f y r(i - 2euv ) ' (22) 

Before we combine these integrals we comment on the identification of UV and IR poles. 
The integrals in Ii and I 2 can be calculated by combining the propagators using Feynman 
parameters and completing the square in the standard manner. I\ is clearly UV finite and IR 
divergent by power counting, and I 2 is clearly UV divergent and IR finite by power counting. 
In I3, if we rescale the loop momenta homogeneously we conclude that the integral is UV 
divergent and IR finite. One may worry if there are IR divergences in J 3 that come from 
k + — > or k~ — > with other components held fixed. In the Appendix, we describe a careful 
evaluation of the integral using contour integration that confirms that the poles are UV. 

It is interesting that the scale Q 2 appears naturally in the evaluation of the integrals Jj. 
The integrals I n and I n and I s are scaleless and cannot know about the scale Q unless it is 
put in by hand. We will see this when we attempt to evaluate these integrals individually 
later in this section. However, after rearranging the integrands, / is expressed in terms of 
integrals, /«, i = 1,2, 3, in which the scale Q naturally appears. It is also important to 
notice that the individual contributions, I i: are free from mixed UV/IR poles thus there is 
no ambiguity in interpreting the e in (fj, 2 /Q 2 ) £ . 

The final result for I is 




21n[^) - I \ / 2m! -4 1 ! 



%v ^uv / \ £ir £m 



(23) 



We must also include the factor of \/~Z^ for each external leg, where 



Z 2 = l-£-C F (± -— ) • (24) 
4vr \e vy em/ 

The final result for the electromagnetic current in SCET to one loop is 



<?(P2)i^ // k(Pi)> = y 1 



, 9 2\n(^-)-3 o 2\n(^)- 
1 + ^C F [4 ^ ±- + ^ J 

4?r \ £ vv £ uv £ir £ir 



.(25) 
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The UV poles are canceled by counterterm for the SCET effective current, 



c.t. 



47r' 



3 2 , 
1 In 



-uv 



£uv £uv 



t< 2 



(26) 



The IR poles are exactly the same as in the full QCD calculation, which is given by 



(q(P2)\f\q(pi)) =7^ 



a 



2 In 



1 + -t^Cf -^r 



-3 



4tt 



^IR 



InM +31n -V 



7T 2 

- 8 + y 

(27) 

and the matching coefficient of the full QCD current onto the SCET one is just the finite 
part of Eq. (T5j 



C{Q 2 /^ 



1 



-\n A +31n 



n 2 ' 



(28) 



from which the anomalous dimension of the effective current 71 is obtained: 

dlnC 



<iln/i 



7i 



with 



7i 



-41n(f 
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(29) 



(30) 



The above results for C(Q 2 / fi 2 ) and 71 where first obtained in Ref. [161 ] . where offshellness 
was used to regulate the IR divergences. Next we attempt to directly evaluate the the 
integrals I n and I s . It is possible to unambiguously determine the double poles in e in these 
graphs but the single poles are ambiguous. Consider first the naive collinear contribution 
I n , which can be written as 



-2ig 2 s C F (fi 2 y 
d d k 



d d k 



V 



+ 



(27r) d [-2pk- + k 2 + iQ][-k+ + i0}[k 2 + z0] 
I 



(2vr) d 



-2pk~ + k 2 + iO][k 2 + i0] 



(31) 



To calculate the first integral in Eq. (EH), begin by performing the k + integral using contour 
integration, then do the d — 2 dimensional integral over the transverse momentum. One is 
left with an integral over k~ which is proportional to 



dk-{2pk~ 



-l-e 



Q 



-2c 



^UV ^UV 



(32) 



In this formula, we have rescaled k~ —>■ pk~ then used a standard result in dimensional 
regularization. The rescaling is required so that that the equation is dimensionally correct 
in d dimensions but the dimensional quantity that is raised to — 2e power is arbitrary and 
we have put the scale Q into the integral by hand. The second integral in Eq. ( 13TI) can be 
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easily calculated by combining the integrals using standard Feynman parameterization and 
the result is 



2r 



d d k 



(2n) d k A 16tt 5 



Q 2 



1 

£lR 



(33) 



where again the scale Q has been inserted to make the result sensible on dimensional grounds. 
The final result is 



a. 
Air' 



■Gi 



Q 2 



2 2 \ 



2 

Sir 



(34) 



There are two sources of ambiguity in the evaluation of this integral: the scale Q has been 
inserted by hand, as discussed earlier, and furthermore, the expansion in e is ambiguous 
because of the mixed l/(euv^m) pole. Because of this, our result for I n (and I fl ) should 
really be regarded as a prescription for defining the integral. A similar situation arises in 
the evaluation of I s , as we will see below. 

For the soft diagram, the contour integration fails and one has to perform the integral 
differently. Inserting by hand a scale in the integral, we rewrite I s in the form 



d d k 



(2n) d (p 1 -k)(p 2 -k)k 2 ' 



(35) 



where, as before, p\ (p 2 ) stands for the momentum of the incoming (outgoing) parton with 
Piipi) = S2 an d Pi = 0- Using the following identity 



1 



1 + 



k 2 



2pi ■ k 



the integral becomes 



with 



Pi ■ k (pi + k) 2 
I s = -2tg 2 C F (fi 2 ) £ Q 2 x [/ S)1 



1,2 



L.2 + Is 



d d k 



1 



{27r) d (px + k) 2 {p 2 + k) 2 k 2 ' 



(36) 



(37) 



(38) 



d d k 



1 



and 



{2n) d { Pl + k) 2 {p 2 + k) 2 
d d k k 



1 1 

+ 



2pi ■ k p 2 ■ k 



(2n) d (pj + k) 2 (2 Pl ■ k)(p 2 + k) 2 (p 2 ■ k) 



(39) 



(40) 



The calculation of I Sj % is straightforward. Multiplying the result with the pre-factor in 
Eq. (|37|) and denoting the result J Sj i one obtains 



a 



47T 



2 \ £ 



P 2 " + "fi" 
£ IR . 



(41) 
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For I Sj 2 and J S: 3, the first step is to combine 

1 1 1 



using the identity 



2 Pi -k(pi + k) 2 2 Pi ■ k k 2 + 2 Pi ■ k ' 



1 f°° 1 

b(a + b) = ] 1 dX [a + Xb} 2 ' 



(42) 



After this step the evaluation of the integral is straightforward and we find 



j _ a s r ( V 2 x 

/s < 2 " T^ Cf {-q 2 



7T 



^UV^IR 3 



and 



1 8,3 





' 2 7T 2 " 




. 4v 6 . 



(43) 



(44) 



The combined result gives 



1 



(45) 



Again, J s is ambiguous because the scale that compensates \i is arbitrary and because of 
mixed 1/(£uv£ir) poles. When our results for I n , I n and I s are combined according to 
Eq. ( 1T81) . the mixed 1/(£uv £ ir) poles in I cancel. We can then separate the UV and IR 
divergent terms, expand (fi 2 /Q 2 ) e in each term and recover Eq. (12"3"|) . The prescriptions 
for defining I n and I n and I s are thus justified a postieri by the requirement that SCET 
reproduce the IR divergences of QCD. 

Our main result for this section is that the soft subtraction in Eq. ([3]) gives the same 
result as the zero-bin subtraction when the one-loop graphs for the jet and soft functions 
are evaluated using DR to regulate the IR as well as UV. While the one-loop integrals in 
the evaluation of the collinear and soft functions are ill-defined, it is possible to see the 
equivalence of zero-bin and soft subtraction at the level of the integrands. The sum of one 
loop collinear and soft graphs is well-defined and reproduces the IR divergences of QCD. 
We gave a prescription for evaluating these integrals which reproduces these results. 



III. DIS AS x 1: SOFT SUBTRACTION 

In this section we extend the analysis of the quark form factor and zero-bin subtractions 
presented in Section II to the DIS non-singlet structure function F 2 (x, Q 2 ) in the threshold 
region, x — > 1. We follow the notation of Ref. [19( with x = Q 2 /2pi-q and p\ is the momentum 
of the incoming parton. We define all quantities to be gauge invariant. All Wilson lines are 
defined on the light-cone and as before we regularize both IR and UV divergences in pure DR. 
It will be shown at 0(a s ) that the zero-bin contributions exist in Feynman diagrams with 
real gluon emission that contribute to the naive collinear matrix elements. For DIS these 
are the well-known jet function and PDF to be defined below. Moreover we will see that 
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the zero-bin is equivalent to the soft contribution in pure DR. When eliminating the double 
counting by generalizing Eq. ([121 to take into account a product of two electromagnetic 
currents we recover, in the SCET formalism, the factorization theorem given in Eq. ([3]) that 
holds to all orders in perturbation theory. 

The general strategy we advocate to eliminate double counting in SCET factorization 
formalism for a given physical quantity is to first decouple the soft gluons from collinear 
SCET Lagrangian by performing field redefinition: £ — > and — > A^'^. In what 
follows, we drop the superscript (0) with the understanding that decoupling the soft gluons 
is already performed. Then one defines naive collinear matrix elements and includes in the 
SCET Lagrangian the zero-bin fields. This enables us to perform perturbative calculations 
extended over all momentum space which include contributions from the soft momentum 
region. These contributions should then be eliminated by performing the zero-bin subtrac- 
tions. We will show at one-loop order that these zero-bin subtractions are equivalent to 
dividing by matrix elements of soft Wilson lines, confirming Eq. The validity of this 
procedure can be justified to all orders in perturbation theory and to lowest order in A 
following the arguments of Ref. 



Let us start by defining the soft factor S which is given by, 

p f dX 



S(l-x) 



i\(l-x)p, 



N r / 2vr 



(0\Tr[Y n (Xn)Yl(Xn) x Fj(0)Y n (0)]|0) 



where p = Q/y/2. The pre-factor is chosen to normalize the leading contribution to 5(1 — x). 
The Feynman diagrams with real gluon emission that contribute to S(l — x) are given in 
Fig. [2J Figs. [2(b) and [2(c) are identically zero due to n 2 = n 2 = respectively. The 
contribution from Fig. |2j(a) is 

S^(l-x) = 2a s C F (fi 2 ypj J**_ 9 .J K J (k 2 )6(k + -(l-x)p) 



(27r) d - 2 k+k- 



6(l-x)+D (x) 



£uv 



-) 



— . (46) 



and the Q 2 -dependence is again put by hand as the integral is scaleless. We use Di[x) for 
the "plus" distributions: 



DAx) 



ln l (l -x) 
1 — x 



0,1. 



(47) 



In going from the first to second line in Eq. fj46|) we have carried out the integral over k + 
first, then integrated over the transverse momentum, and finally integrated k~ from and 
oo. Taking the contribution from the mirror diagram of (a) and adding 21 s x 5(1 — x) (I s 
is given in Eq. (J3SD ) to include the virtual contributions we get the soft factor to 0(a s ) in 
pure DR: 



5(1 -x) = 5(1 -x) 

+ ^C F 
Air 



(48) 



Q 2 



^) x4 



+ 



-uv 



^UV^IR 



5(1 -x) + D (x) 



£uv 



1 



The renormalized soft function to 0(a s ) is 

S R (l-x) = 5(l 



x 



-4) 



Do(x) 
£ir 



(49) 
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(a) (b) (c) 

FIG. 2: Real gluon contribution to the soft factor. 



Now we consider the PDF. Our definition of the PDF is analogous to the standard one 
in QCD j20[] however it is expressed in terms of the naive collinear SCET fields, 

m = \ J ^e^(P||jAn)iy n (oo,0;An)x 7 + ^t (TO)0 . 0) ^ (0) | P) . (50) 

The Feynman diagrams with real gluon emission that contribute to the partonic PDF are 
given in Fig. [31 The contribution from Fig. [3](c) is again identically zero due to n 2 = 0. The 




FIG. 3: Real gluon contribution to the PDF. 



contribution from Fig. E](b) is nonsingular in the limit x — > 1. Its contribution should be 
omitted since it is subleading in the SCET expansion parameter, A, which for DIS in the 
threshold region has to be taken as A 2 = 1 — x. The remaining contribution from Fig. [3]^a) 
is given by: 



= 2a s C F {^y J { ^yJ k X 5{k2)5{k+ ~ (1 " X)P) 

2 2 \ 



2 \ e 



-— 5{1 - x) + D {x) 



(51) 



Comparing the last result with the contribution from Fig. E^a) given in Eq. ( H6l) we see that 
in the x — > 1 limit both contributions are identical. 

For the PDF the leading contribution in the x — > 1 limit (or equivalently in A) comes 
entirely from the zero-bin region as we can see from the second 5-function in Eq. (I3T1) . k + 
which is supposed to be 0(Q) (since k is collinear to the incoming parton momentum with 
Pi = p = Q/y2) is restricted to be equal to (1 — x)p which is 0(Q\ 2 ). Since k~ is also 
0(QX 2 ) then the first 5-function enforces the transverse components to be 0(QX 2 ). Thus 
the contribution comes from the zero-bin (soft) region. Notice that Fig. E^a) is the only 
diagram relevant to the PDF as x — > 1 because we have defined our PDF with the naive 
collinear fields, which allowed us to integrate over all momentum space including the zero-bin 
region. If we had restricted the SCET definition of the PDF to contain the purely collinear 
fields (i.e., so k + is 0(Q)) then Fig. [3](a) would not contribute as the second 5-function in 
Eq. fl5Tl) could not be satisfied by power counting arguments. Existing SCET treatments of 
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DIS in the x — > 1 limit differ in their treatment of Fig. [3](a) [17|, [18|] . Ref . [17| drops this 
diagram and argues that collinear modes cannot contribute because of kinematic constraints. 
This is essentially equivalent to our argument above that the diagram receives support only 
from modes whose scaling is soft rather than collinear. In our approach this contribution 
is removed from the purely collinear jet function by a zero-bin subtraction. Our purely 
collinear PDF (i.e. the naively evaluated PDF minus the zero-bin subtraction) is analogous 
to the function gp that appears in the factorization theorem of Ref. [18|. However, zero- 
bin subtractions are required for both the PDF and the jet function, as we will see below, 
whereas Ref. 18[ only includes a zero-bin subtraction in their evaluation of the PDF. 

The final state jet function represents outgoing (in the n direction) collinear partons with 
invariant mass Q 2 (l — x), which is assumed to be much larger than Aq CD . To define the jet 
function, we start with the following two-point correlation function of collinear fields [9| 



(0\T M(z)Uz)U0Wn(0) \0) = i-t= 



d 4 k 



-ikz 



V2J (2tt) 



J(k) . 



(52) 



The naive dimensionless jet function, to be denoted by J(Q 2 , x), is related to the absorptive 
part of J and normalized at leading order to 5(1 — x), so that J(Q 2 ,x) = x ImJ'. 
This definition is the gauge invariant version of the jet function defined in Ref. [3( with the 
full QCD fields replaced by the naive SCET ones. The Feynman diagrams with real gluon 
emission that contribute to the jet function are given in Fig. HI 




(a) (b) (c) 

FIG. 4: Real gluon contribution to the outgoing jet function. 



The contribution of Fig. [J](c) is zero due to n 2 = 0. From Fig. H^a) we have 

d d k 8pl~ 



J( a \Q 2 ,x) = a s C F (ii 2 Y 



(2tt; 



d-2 



k-V 



— k~)5(k )5((l — k) ) 



(53) 



Q_ 

v/2 



where / is the momentum for the outgoing jet, whose components are Z + = -^(1— x); I' 
and l± = 0. k is collinear to the outgoing quark momentum with k~ = Q, k + = QX 2 and 
k± = QX. Recalling that in threshold region for DIS we identify 1 —x = X 2 then simple power 
counting shows that the above contribution scales as 1/A 2 . Carrying out the integrations 
over k + , k±, and k~ (in that order) we find 



J {a \Q\x) 



Air 



C F x 4 




5(1 -x) - 



Do(x) 



-D (x) + D 1 (x) 



\\ 5(1-*)+ m \%)Do(x) 



.2 \^ 

The contribution from Fig. H^c) can be obtained in similar manner and is given by 



J {c) (Q\x) = ^C F 



5(1 -x) 

^IR 



l-\n[^))5(l-x) 



Do(x) 



(54) 



(55) 
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Including the contribution from the mirror of Fig. H](a) we get for the real gluon contribution 
to the jet function 



>^real(Q > 




X 



AD (x 

SIR 



+ 4 J D 1 (x) 



(56) 



- 3 In 



+ 7 



7T 



5(1 



Let us now consider the zero- bin contribution included in the above result. For Fig. HJa) we 
take the gluon momentum k to be soft. Then we ignore k~ relative l~ in the numerator of 
Eq. ( 1531) and we drop k 2 in the last 5-function as it scales as A 4 compared with I 2 = Q 2 X 2 . 
We then get a contribution that also scales as 1/A 2 like in the collinear region: 



J^(Q 2 ,x) = 2a s C F (fi 



2\e 



Tn° F {~Q 2 



d d k 



-5(k 2 )5(k + - V 



(2tt)^ 2 k 

-—8(1 -x) + D (x) 
Sir 



Suv 



2 

Sir 



— , (57) 



The last result is exactly the same as the ones given in Eq. (1461) and Eq. (!5"Ti) (when taken 
to the x — > 1 limit.) Noticing that the zero-bin contribution from Fig. H](b) is subleading in 
A so the soft contribution to the real gluon emission for the jet function comes only from 
Fig. Ha). 

We now include the virtual contributions to the soft factor, PDF and jet function. For 
the soft factor the only contribution comes from Fig. [0(c). The sum of the real and virtual 
contributions is 



S(l 



5(1 -x) + —C F x 4 

47T 



suv 



-i-) 

sir/ 



suv 



-5(1 -x) + D (x) 



(58) 



Note that the soft function is scaleless and the final answer is proportional to 



j i_ 

The factor (fi 2 /Q 2 ) 6 can be set to unity since when we expand this factor in powers of e 
and mulitply by — the finite logarithms cancel. This makes physical sense since this 

quantity is scaleless. A similar situation arises for the PDF, 4>(x), however, the jet function 
will have logarithms of Q because it knows about the scale Q\/l — x. The renormalized soft 
function to 0(a.) is 



S R (1 



X 



5(1 



x) + —C F x (-4) — — 



47T 



Sir 



The large N moments of the soft function are 



(59) 



4:71 



Sn = 1 + x 4 



In AT 

Sir 



(60) 



with InN = Ne 113 . For the PDF, after including virtual contributions, we find 



4>(x) 



6(l-x) + -- t C F 

47T 



SUV 



Sir/ 



8(1 - x) + AD (x) 



Suv 



1 

Sir 



(61) 
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The renormalized PDF is 



4> R (x) = 5(1 -x) + ^C F ( - 



and in moment space we get 



Arc 



^IR 
1 



[35(1 -x)+ 4:D (x)} 



[3-41nJV] 



(62) 



(63) 



which is a well-known result. The anomalous dimension of the (naive) PDF can be imme- 
diately read off from the UV poles in Eq. flBTl) : 



a 



l2,N 



—C F x 2[41niV-3l . 

47T 



(64) 



The virtual contribution to the jet function comes from Fig. [0(b), its mirror diagram, 
and the wave function renormalization. In the result for I n (which is equal to I n given in 
Eq. (|34"|) ) we expand the factor (fj, 2 /Q 2 ) e and get for the naive jet function 



J(Q 2 ,x) 



o:. s 



An 



AD (x) 



+ 2 In 2 



Q 2 



- 3 In 

4 



^ir 

/i 2 



4 J D 1 (x)+41n(^-)D (x) 



+ 7 



7T 



5(1 -x) 



SUV^IR SUV 



5(1 -x) . 



(65) 



Next we have to include the zero-bin subtraction which as we have seen earlier is equivalent 
to subtracting the one-loop soft contribution. It is easy to see the effect of the zero-bin 
subtraction is to replace 1/erR everywhere in Eq. (1651) with 1/euv- The remaining UV 
divergences can now be removed by counterterms and the renormalized jet function is 



a. 



J H (Q\x) = 5(l-x) + -^C F 

47T 



AD 1 (x)+A\n[%)D (x) 



+ 2^ "Sin ^ +7-^ 5(1 



Q- 



X 



(66) 



The results presented so far for the zero-bin contributions (real and virtual) for the naive 
PDF and jet function show that the renormalized soft factor has to be subtracted from each 
one of these functions to obtain the truly collinear contributions. Thus we find that the 
renormalized collinear matrix elements are 





which to 0(a s ) are given by 



1 



(67) 



1 + t= C F x (-3)— , 

47T Sir 



(68) 
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and 



Nn 2 ) \Nfi 2 J 3 



(69) 



The last result for the collinear jet is finite and is equal to the matching coefficient for DIS 



at the intermediate scale, fi 2 = Q 2 /N, in the analysis of DIS as x — > 1 in Ref. 16) . Thus 
the factorization theorem for DIS in the threshold region reads 

F^ = H(Qy^)jU N S§ = H(Qy^)(^j f^j S §. (70) 

where if (Q 2 //i 2 ) is the square of the matching coefficient C(Q 2 j /i 2 ) given in Eq. (|28|) . With 
the above results for J§, (fi N and S§ we get 

H(Q*/^)J§<f>*S« = 1 + ^cJ- — [3-41niV]-21n 2 f^+61nf^ 

4vr [ £ IR \H J \H 



N^ 2 J \Nfj, 

which agrees with the moments of DIS structure function in the large- N limit calculated in 
full QCD. It is straightforward to show that the results in Eq. (!67|) which were shown to 
hold to 0(a s ) in pure DR can be obtained by performing field redefinitions [l3[ on the naive 
SCET fields and one obtains the naive collinear matrix elements divided by the soft Wilson 
line matrix elements as in Eq. ([3]). 

Based on the factorized form of the non-singlet DIS structure function given in Eq. (ITUj) 
we now comment on the resummation of the large logarithms in the threshold region for 
DIS. In moment space Eq. (1701) can be written in the following form (henceforth we drop 
the superscript R with the understanding that we consider only renormalized quantities) 

F 2 ,n = H{Q 2 /f)J N (J^ <M/i 2 ) . (72) 

The hard part H depends only on Q 2 //i 2 and it is obtained by matching the full QCD 
current onto the SCET one (at the higher scale Q 2 ) order by order in perturbation theory. 
The anomalous dimension of the SCET current 71 is then calculated from the matching 
coefficient (through Eq. (I2"9"j) ) and is used to run down to the intermediate scale of DIS: Q 2 / N. 
The quantity J^, which is IR safe, depends only on the intermediate scale as can be seen in 
Eq. (l6"91 . Thus the first two terms on the right hand side of Eq. (IT2"j) are perturbative and 
IR safe. Below the intermediate scale we are left with only one non-perturbative quantity, 
the PDF taken to the large N limt. By exploiting the standard Altarelli-Parisi kernels 
with anomalous dimension 72 (taken to the large iV-limit) we can then evolve the PDF to 
an arbitrary factorization scale. The two stage running between Q 2 and Q 2 /N with 271 
and between Q 2 /N and some arbitrary factorization scale //p with 72 resums all the large 
logarithms in moment space. This has been established to be equivalent to the standard 
pQCD resummation [I, [2l[ in Ref. [22| to all orders in the strong coupling constant and to 
arbitrary sub-leading logarithms. 
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IV. ZERO-BIN AT HIGHER ORDERS 



In this section we consider the abelian virtual two-loop diagrams that contribute to the 
n-collinear jet function (0\Wn°^ ffi \q(p)} that appears in the factorization theorem for the 
quark form factor. We will show how dividing by the soft factor reproduces the zero-bin 
subtraction for these diagrams. 

The abelian two-loop diagrams are shown in Fig. Before considering these in detail 
we make some general comments about the zero-bin subtraction at two loops. Let k\ and 
k 2 be the loop momenta, which are routed so that k\ and k 2 correspond to the virtual 
gluon momenta, since the zero-bin associated with any fermion lines is easily checked to be 
subleading in the A expansion. Since we integrate over all momentum space, we need to 
subtract the contribution where either of the ki scales like a soft momentum rather than 
collinear. In Fig. \5[ the region where k\ (k 2 ) is collinear is separated from the region where 
k\ (k 2 ) is soft by the vertical (horizontal) dotted line. The collinear matrix element gets 



ko 



Region-I 



Region-Ill — ► 




FIG. 5: Zero-Bin Regions At Two-Loops. 

contributions only from Region-I. Region-II corresponds to k 2 soft and k\ collinear, Region- 
Ill corresponds to ki soft and k 2 collinear, and Region-IV corresponds to both k\ and k 2 
soft. 

Denote the integrand of the two-loop integral 1{k\ 1 k'^) 1 where the superscript c in ki 
means that ki is assumed to have collinear scaling. The naive two-loop integral is 

( " 2)2 7 S? I ^w,® = J^m,k$. (73) 

The integrands for the zero-bin subtractions for Region-II and Region-Ill are — X(fcJ, k^) and 
—T{k{, k 2 ), respectively, where the superscript s in kf means that the integrand is evaluated 
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assuming that k{ satisfies the soft scaling and the integrand is expanded to lowest order in 
A. The zero-bin subtraction from Region-IV is subtle because this Region has been doubly 
counted both in the original naive collinear integrals and also in the zero-bin subtraction for 
Region-II and Region-Ill. From the original integral we need to subtract a contribution in 
which both k\ and k 2 are s °ft- We will denote the integrand for this zero-bin as X^ikl, 
where L\ denotes the limit when k\ and k 2 are taken to be soft simultaneously. From 
the zero-bin subtraction for Region-II, we must perform a second subtraction that comes 
from kl becoming soft after having first made the soft approximation for k 2 . We call the 
integrand for the zero-bin subtraction defined by this limit lL 2 (kf, /cf). Likewise, we have 
a similar subtraction from the zero-bin of Region-Ill, denoted Xi 3 (kl, k 2 ). In general, the 
integrand Zz^&f , k 2 ) depends on the order in which we take the soft limits, so there are really 
five different zero-bin subtractions in the two-loop calculation. The result for the collinear 
contribution to each two-loop diagram is 



k\,k2 



li^k^ , k' 



T(k 1 ,k 2 ) Xi 2 (k 1 ,k 2 ) 



T(k l ,k 2 ) Xl 3 (ki, k 2 ) 



l Ll (klk s 2 )\ (74) 



The various order of limits defining the integrands Ii i (kl, k 2 ) are shown in Fig. [51 




(a) (b) (c) (d) 




(e) 

FIG. 6: All the two-loop abelian diagrams that contribute the n-collinear jet function. 

Next we turn to the evaluation of the individual diagrams. Figs. [6](a) [6](c) and [6]^d) 
have color factors proportional to C\ while the color factors for Figs. [6]^b) and [6|e) are 
Cf{Cf — Cji/2). The results for an abelian theory are obtained by taking Cp — > 1 and 
Ca — * 0, in which case the color factor is one for all graphs in Fig. [6j 
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Consider first the two-loop SCET diagrams in Fig. [6](a) and[6](b). Working in Feynman 
gauge, using pure DR gulator, and ignoring the ze's in the propagators as they are 

not relevant to our discussion, we get from Fig. E^a), 



JW = -2gi 



(p+ _ k + _ fe+ 



fcl,fc2 



^ (p ^ kl „ k2) 2 k 2 k 2 {kt + ^ 



( P + -kt) + (p" 1 



(75) 



where fc* are both collinear to the incoming parton with momentum p. From Fig. [6](b) we 
have 



i {b) = -2gi 



(p+ - A;+ 



*2 + ) 



fc x ,fca (p - An - k 2 ) 2 k\kl{kt + fc 2 



(p+ - (p-» 



A n 



(p — ki) 2 k 2 (p — k2) 2 kf 



(76) 



In Eqs. (1751) and (1761) . we have chosen to write the integrand so it is symmetric under 
interchange of k\ and ki. 

The integrands for Figs. E^a) and [6(b) in Region- IV are 



z£?(*j,*s) = -^ 



i 



i 



2 k-^k^i^k-^ ~\~ k<2 ) — I - ^2 ) ^i ^i 



i i 

+ 



fc 2 fc 2 



and 



^Li (^1) ^2) — ~dtCp x - 



2 /c-^ i^k-^ ~\~ k<2 1 f^i — 1~~ ^2 



1 2 2 1 



(77) 



(78) 



Note that the contributions in Eq. ( 1771) and Eq. ( 1781) are not subleading in SCET power 
counting. For the two loop momenta in the soft region, the measure of the loop integral 
scales as A 16 and each integrand scales as 1/A 16 thus the contribution is 0(1). The sum of 
these two zero-bin contributions is 



Jk 1: k 2 v 



1 



1 



1 



i,fc 2 2 (k 1 k 1 k 1 )(k 2 k 2 k 2 



(79) 



where I s is the one-loop soft contribution given in Eq. (fT6j) . The two-loop result is consistent 
with the exponentiation theorem for the abelian amplitudes of soft gluon radiation in a web 
(i.e., two soft Wilson lines). 

Now we consider the case when one collinear momentum becomes soft while the other is 
kept collinear. In Eq. (1771) let us take k\ to the zero-bin region. The integrand is then 



2< 8 >(*f,*§)=-2^ 



(P H 



[(p - k 2 ) 2 - 2ki(p+ - k+)}k 2 k 2 k + 



P 



+ 



1 



( P -k 2 ) 2 k+ -2fcffc; 



(80) 



The measure in the loop integral in this case scales as A 12 because there is one collinear and 
one soft loop momentum. The factor outside the brackets in the integrand scales as A~ 8 , 
the first factor inside square brackets scales as A -2 and the second scales as A -4 . Therefore, 
only the second term contributes so we get 



TM(k' 1 ,kZ) = -2gi 



(p+ - ^ 2 



[(p 



k 2 ) 2 -2fcf(p+ -k+)}k 2 k 2 k+ 



1 



2/c ^ k-^ 



(81) 
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Similarly for k\ collinear and k% soft we find 



l^(k c 1 ,k s 2 ) = -2g 4 s 



kt 



[(p-h) 2 

The analogous zero-bins for Fig. [6](b) are 

(p+ - fc. 



2/c 2 (p + &i )]&f 



X( b )(^,^) = -2(7. 



and 



l^ik^k 



(p+ - fc+) 



P 



5) 2 ^[(p-fc 1 ) 2 -2fc 2 "(p 



«, 1 jj^A^fia 



(p - k 2 ) 2 kf _ 

p + — fc^" 
(p — k\) 2 k 2 



(82) 



(83) 



(84) 



Upon summing the zero-bins for each Region, the integrand factorizes. For Region-II, we 

get 



x^{ki,k s 2 )+x^\ki,k s 2 ) = g. 



(p+ - kt) 



X 



(85) 



(r, — hA 2 h 2 h + h 2 h~h + 
yp tvi) hv 1 hv 1 ^2^2 ^2 

and for Region-Ill we obtain the same with k\ and k-i interchanged. The zero-bins for 
Region-II and Region-III combine to give 



fci,fc 2 



(!«(*;, k c 2 ) +i^(ki, k s 2 )+^ b \kt, k c 2 )+i^(k s 2 , k{)) = i n ■ i s 



(86) 



where I n is the one-loop contribution to the n-collinear jet in Eq. ffT4l) . Finally, we have to 
perform the second zero-bin subtraction from the zero-bins corresponding to Region-II and 
Region-III, i.e. the terms with integrands X L2 (kf, k 2 ) and T Lz {k\, k 2 ) in Eq. (173]) . It is clear 
from Eq. (1851) that this is simply |J 2 for each Region. Thus the final result for Figs. E(a) 
and E](b) including all zero-bin subtractions is 



J« + /<»> - [I n ■ I s ~ II 



(87) 



where 1^ and correspond to the naive evaluation of Figs. El^a) and[B](b), respectively. 

Next we turn to the evaluation of Figs. Elc)-(e). Denote the momentum that flows into 
the gluon attached to the collinear Wilson line as k 2 while the momentum flowing through 
the other gluon as k\. For these three diagrams it is easy to show that the only zero-bin 
contribution which is leading in A comes from the region where k\ is collinear and ki is soft. 
In Figs. [6]^c) and (d), k\ flows through the self-energy subgraph. The self-energy in the 
collinear sector of SCET is evaluated in Ref. [9fl with the result that 



a s q 2 jL T(e)T(l-e)T(2-e) 
47r2g+2 T(2 - 2e) 



-q 



where q 2 is the virtuality, which in our case is q 2 = (p — A^) 2 = —2p + k 2 + 0(X 4 
this result into the two-loop graphs in Figs. M^c) and (d) yields the integral 



fci ,fc 2 



(88) 
Inserting 

(89) 



te; r( t )r(i_- t) r(2-o t 



T(2 - 2e) 



j. 2 k 2 k 2 k\ 



(2 P + k-y 
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To evaluate Fig. [6](e), we perform the k\ integral which is simply the collinear vertex correc- 
tion in the limit that the external gluon is taken to be soft. Equivalently, this subgraph is 
simply the one-loop vertex correction to the coupling between the collinear quark and soft 
gluon. The result after doing this integral is 

t i^(k c k s ) - , a , 2 r(g)r(i-6)r(2-6) 4 r 1 
Lto ( " 2) " s r(2 - 2e) (e " } J k2 ktk~ 2 kl {2p k > } 

We find that the zero- bins for Fig. E](c),(d) and (e) add up to zero. The cancellation can be 
partly understood on the basis of the QED Ward identity which relates the UV divergent 
pieces in the self-energy and vertex correction. 

There are no zero-bin subtractions beyond those given in Eq. (|87|) . We now compare this 
result with the soft subtraction, Eq. (fT2j) . which yields 

(0\Wlin\q(pi)) = l + In + I {a) ~ (e) + ... (QQ) 

(0\YlY n \0) " l + / s + i/ s 2 + ... 

= 1 + (I n ~ Is) + (/ (a) - (e) - In ■ Is + \l 2 )j + - • 

In the right hand side of the first line, we have expanded the numerator and denominator to 
O(a^), separately. The one-loop naive collinear integral is /„, j( a )~( e ) is the naive collinear 
evaluation of the sum of graphs in Fig. El I s is the one-loop soft integral, and the second 
order term in the denominator is a consequence of the exponentiation theorem of the abelian 
contributions inside a web mentioned earlier. In the second line, we have expanded the quo- 
tient to second order in a s . The 0(a 2 s ) term is consistent with Eq. (IH71) and the cancellation 
of the zero-bin in the sum of Figs]6^c),(d) and (e). Therefore, dividing by the soft Wilson 
line and the zero-bin subtraction give the same result. It would be interesting to extend the 
analysis of this section to the non-abelian theory. 



V. CONCLUSIONS 

In this paper we have considered the two prescriptions deivsed to remove overlapping 
contributions to collinear matrix elements, the soft and zero-bin subtractions. We have 
demonstrated explicitly the equivalence of the two prescriptions for the abelian contributions 
to the quark form factor up to two-loop level using DR to regularize both UV and IR 
divergences. We also studied DIS in the threshold region to 0(a s ) in SCET. In our treatment 
for DIS all Wilson lines were defined on the light- cone. The essential result is that soft 
contributions to naively defined collinear matrix elements have to be subtracted in order to 
derive proper factorization theorems as in Eq. ([3]). We have shown by explicit calculation 
that soft and zero-bin subtractions are equivalent in the examples studied in this paper. 
Our results obtained by fixed order pertubative calculations can be extended to all orders in 



the strong coupling by performing field redefinitions as proposed in Ref. [13j| with a suitable 
choice of IR regulators. 
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VI. APPENDIX 

Here we explain in detail why the poles in Eq. ([22"]) are UV using contour integration. 
Let us first integrate over the light-cone component k~ with contour integration. The poles 
are 

,-_ \k±\ 2 iO ,__ \k ± \ 2 + 2pk+ *0 

k ~ 2(k+-p) 2(k+-py * - 2^ 2F' * "' U - W 

There are two regions that contribute, < k + < p and k + > p. For < k + < p we have 
to pick the second pole. For k + > p we have to pick the third pole. After performing 
the contour integration we integrate over \k±\ 2 in d = 2 — 2e. This will introduce the first 
r(exrv)- We then integrate over k + . For the contribution from < k + < p, the integration 
over k + will introduce another pole as k + — > 0, r(— e). This pole has to be taken as UV 
for the following reason: When k + — > 0, then from the second pole in Eq. (1911) we see that 
k~ — > oo. This combination clearly means that both k° and k 3 approach infinity (in addition 
to the UV pole from \k±\ 2 — > oo). The same reasoning applies for the contribution from 
the region k + > p where we get r(— e) from k + — > oo. The third pole in Eq. (I9TI) indicates 
k~ — > so we again have one light-cone component approaching zero while the other one 
approaches infinity. Thus both contributions have a double UV poles and the total result 
is given in Eq. (I2"2"j) . The observation that a vanishing light-cone component may lead to a 
UV divergence (as opposed to IR divergence) was also discussed in Ref. [141 ]. 
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